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We report new constraints on flavor-changing non-standard neutrino interactions (NSI) using data 
from the MINOS experiment. We analyzed a combined set of beam neutrino and antineutrino data, 
and found no evidence for deviations from standard neutrino mixing. The observed energy spectra 
constrain the NSI parameter to the range —0.20 < e Mr < 0.07 (90% C.L.). 
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It is well established from solar, atmospheric, reactor 
and accelerator experiments [U-Q that neutrinos undergo 
flavor change as they propagate. This phenomenon can 
be explained by the quantum mechanical mixing of neu- 
trino flavor and mass eigenstatcs. The mixing can be 
parametrized by three angles, #12, #13, #23, and a CP- 
violating phase, S Q. The standard neutrino oscillation 
mechanism requires that at least two of the three ac- 
tive neutrinos are massive. While the phenomenon of 
oscillations can occur in vacuum, the presence of mat- 
ter along the neutrino path allows for alternative fla- 
vor changing mechanisms such as the Mikheyev-Smirnov- 
Wolfenstein (MSW) matter effect @. This effect alters 
the survival probability of electron neutrinos propagat- 
ing through matter since electron neutrinos can have ad- 
ditional interactions with electrons in the surrounding 
medium, but does not affect the survival probability of 
muon or tau neutrinos. 

Non-standard interactions (NSI) [l(J could occur be- 
tween muon or tau neutrinos and matter, and could al- 
ter the flavor content of a neutrino beam as it propagates 
through the Earth's crust in a manner similar to standard 
matter effects. Searches for NSI have already been per- 
formed with atmospheric neutrinos (llj . However, non- 
standard matter effects are, in general, different for neu- 
trinos and antineutrinos. Accelerator-based oscillation 
experiments offer a powerful tool to search for NSI with 
their ability to produce well-understood beams of neu- 
trinos and antineutrinos separately. Furthermore, MSW 
and NSI effects depend on the neutrino's path length in 
matter, and beam neutrinos travel a well defined distance 
to the detector. Recent papers have discussed the com- 
patibility of NSI with MINOS data [HE1. This Letter 
describes the first direct search for NSI which simulta- 
neously fits the separate energy spectra of neutrinos and 
antineutrinos in a long-baseline experiment. 

Short-baseline neutrino experiments have explored and 
constrained non-standard interactions fl4| . In such 
experiments both charged-current (CC) and neutral- 
current (NC) NSI can be studied. While long-baseline 
experiments can also constrain interactions with matter 
directly in either near or far detectors, their improved 
sensitivity arises from using the Earth's matter along the 
long neutrino path as the interaction medium. Charged- 
current NSI produce a final state chargcd-lcpton which is 
absorbed by the Earth's matter and not observed. How- 
ever, NC NSI that occur along the neutrino path will 
produce a final state neutrino of a different flavor, alter- 
ing the flavor content of the beam [l3| • We consider only 
NC NSI in this study. 

The NSI Hamiltonian can be included as a perturba- 
tion to standard oscillations. This Hamiltonian is propor- 
tional to the matter potential V = \Z2GfN s , with Gf, 
the Fermi coupling constant, and 7V e , the electron density 
in matter, analogous to the MSW matter effect 0]. In 
general, the Hamiltonian has both flavor-conserving and 



flavor-changing components. In flavor-conserving NSI, 
the NC scattering between the neutrino and matter does 
not alter the neutrino flavor. Flavor-changing NSI, on 
the other hand, do not conserve lepton flavor number; 
the final state neutrino is in a different flavor eigenstate 
from the initial neutrino. In a disappearance experi- 
ment, flavor-changing NSI have a greater effect on the 
flavor transition probability than flavor-conserving NSI 
due to interference between amplitudes of standard and 
non-standard matter interactions 151. We consider here 



only flavor-changing NSI and set to zero flavor-conserving 
amplitudes to which MINOS has no sensitivity. 

The resulting NSI Hamiltonian in the two-flavor ap- 
proximation is 



hnsi — v 
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where the coefficient gives the strength of the NSI ef- 
fect on transitions between /x and r flavors. We only con- 
sider the real part of e MT , which is sensitive to differences 
between and F M survival. We define the vacuum oscil- 
lation length for neutrinos of energy E as Lq = ( Am 2 ) ■ 
The difference between the squares of the second and 
third neutrino masses, A? 
the same 
cillations [16 
defined as 
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arameter that governs standard neutrino os- 
The NSI matter oscillation length is 
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[1 ± 2sm(28)L e^ T \V\ + (L £^ M) 2 ] 5 
The survival probability can then be written as 
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with mixing angle 8 = 623, and neutrino path length L. 
Standard oscillation parameters 8 and Am 2 arc taken 
to be the same for neutrinos and antineutrinos. The ± 
signs in Eq. ([2]) arise from the matter potential, V , which 
is positive for neutrinos and negative for antineutrinos. 
The parameter e^ T is real- valued and carries its own sign. 
A positive value of e^r implies that the neutrino disap- 
pearance probability is greater than the antineutrino dis- 
appearance probability, and vice versa. 

The MINOS experiment measures disappearance of 
muon neutrinos and antineutrinos in the NuMI beam [l8[ 
using two detectors. The event energy spectrum of the 
low energy NuMI beam, used in this measurement, peaks 
at approximately 3 GeV [l!|. Its focusing components 
can be tuned to produce a beam with an event composi- 
tion of 91.7% 7% Vp, and 1.3% (v e + v e ) in neutrino 
mode, or of 58% v^, 40% F M , and 2% (v e + V e ) in an- 
tineutrino mode [20j. The Near Detector (ND), with a 
fiducial mass of 23.7 tons, measures the neutrino and 
antineutrino energy spectra 1.04 km downstream of the 
production target. The Far Detector (FD) is located in 



3 



the Soudan Underground Laboratory and has a 4.2 kilo- 
ton fiducial mass. It measures the energy spectra 735 km 
downstream of the production target. Both detectors are 
magnetized steel-scintillator tracking-sampling calorime- 
ters designed to measure the energy and the sign of the 
charge of muons produced by and 77^ interactions. In 
each detector, muon neutrino and antineutrino CC in- 
teractions are separated event-by-event using the sign of 
muon track curvature. 

The results presented here are based on an exposure of 
7.09 x 10 20 protons on target (POT) in neutrino mode, 
combined with a 2.95 x 10 20 POT exposure in antineu- 
trino mode. MINOS previously reported a two-flavor os- 
cillation analysis on this neutrino dataset [Til 17 1 and 



the first direct measurement of antineutrino oscillation 
parameters from the antineutrino sample (20| . Due to 
the opposite sign of the matter potential in Eq. ([3]) for 
neutrinos and antineutrinos, NSI, if present, will alter 
the survival probability of neutrinos and antineutrinos 
in opposite directions. The magnitude of e^ T is propor- 
tional to the difference in probability between neutrinos 
and antineutrinos, and the sign of £ MT is determined by 
the sign of the probability difference. 

We select and CC events inside the fiducial vol- 
ume by identifying interaction vertices with a muon track 
and possible hadronic activity. The neutrino energy is 
reconstructed by summing the muon track and hadronic 
shower energies. Muon energy is measured using range 
for muons that stop in the detector and curvature for 
muons that exit. The hadronic energy is determined us- 
ing a fc-Nearest Neighbor (fcNN) technique 2l|. We re- 
quire the muon charge be negative for and positive for 
F M events. 

To reject NC interactions we use a discriminant ob- 
tained by combining four event characteristics into a fcNN 
variable [22| . The selection criteria optimize selection ef- 
ficiency and sample purity to obtain maximum sensitivity 
to oscillations [2l|, |23|. Far Detector selection efficiencies 
for neutrino and antineutrino samples arc 93% and 97%, 
with purities of 99% and 94%, respectively. Because the 
neutrino sample has a larger component of highly inelas- 
tic events, this optimization process reduces the optimal 
neutrino selection efficiency in favor of lower NC back- 
ground. The lower overall purity of the antineutrino sam- 
ple results from the much larger neutrino contamination 
at higher energies in the antineutrino mode; however, in 
the region of interest to oscillations and NSI the contam- 
ination is smaller. 

The FD neutrino and antineutrino spectra in the ab- 
sence of flavor change are predicted using the ND data by 
first correcting the ND spectra for inefficiency and back- 
grounds and then extrapolating to the FD by a trans- 
fer matrix obtained from simulation 0, [24| . We predict 
2073 neutrino and 273 antineutrino events without oscil- 
lations, and observe 1654 and 193 events, respectively. 

The neutrino and antineutrino energy spectra are fit 
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FIG. 1: Far Detector distributions of selected neutrino (top) 
and antineutrino (bottom) events. Black dots represent data, 
the dashed histogram shows the prediction in the absence of 
oscillations and the solid histogram shows the prediction for 
the model in Eq. ((SJ) at the values obtained in our fit to the 
oscillation and NSI model. 



simultaneously to three parameters, |Am 2 |, sin 2 28, and 
£ MT , in the combined oscillation and NSI model in Eq. 
using a binned log-likelihood. The value of the mix- 
ing angle is constrained to be physical by asserting 
< sin 2 (2f5) < 1. The resulting simulated energy spec- 
tra, obtained by using the best fit values, are shown in 
Fig.[T]supcrimposed on the full neutrino and antineutrino 
spectra. 

The overall systematic uncertainty in the measurement 
is much smaller than the statistical uncertainty. How- 
ever, the difference in the relative sizes of the neutrino 
and antineutrino event samples results in a large differ- 
ence in their statistical uncertainties, while the system- 
atic uncertainties are comparable. Systematic uncertain- 
ties are included in the fit to ensure that the neutrino 
sample does not have a disproportionate impact. The 
largest sources of systematic uncertainty are the muon 
energy scale, the hadronic energy scale, the NC back- 
ground, and the relative normalization between the Near 
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and Far detectors. The muon energy scale uncertainty 
is 2% for events where muon energy is determined from 
range and 3% where the energy is obtained by measuring 
the curvature in the magnetic field. 

The uncertainty on the measurement from range is 
based on a comparison of tabulated muon spectrum 
power in MINOS detectors with tables in [25|, and in- 
cludes our uncertainty in the material composition of 
the detectors. The curvature uncertainty was found 
by comparing the curvature measurement to that from 
range for stopping muons 16j. The hadronic energy 
scale uncertainty is energy-dependent and smaller than 
10% integrated over all hadronic energies. Its multiple 
components arise from uncertainties in hadron produc- 
tion modeling and from the uncertainty in detector re- 
sponse The NC background is less than 2% of 
both the neutrino and antineutrino samples integrated 
across all energies. The size of its uncertainty, which is 
dominated by hadronic shower modeling, is estimated to 
be 20% by comparing experimental data to simulation. 
The 1.6% normalization uncertainty at all energies arises 
from modeling reconstruction differences between Near 
and Far detectors, and the uncertainty in the number of 
protons on target fl(| . 

These four sources of systematic uncertainty are in- 
cluded in the fit using penalty terms. The best fit pa- 
rameters from this procedure are found to be 



Am 2 = 2.39±fH? x 10~ 3 



sin 2 (2fl) =i.oo±8:8g, 



eV 2 



Ear = -0.07 



+0.08 
0.08' 



with the allowed region -0.20 < < 0.07 (90% C.L.). 

The penalty terms from systematic uncertainties have 
a negligible effect on the fit; each penalty term pulls 
the best fit point by much less than one standard de- 
viation. The allowed regions of fit parameters are shown 
in Fig. [5J where three two-dimensional slices from a 3D 
likelihood surface are chosen by marginalizing over the 
third parameter. The obtained oscillation parameter val- 
ues are i n good agreement with previously published re- 
sults in 



20l | . Within errors the fit is consistent with 
no contribution to flavor change from NSI. This conclu- 
sion is in agreement with recent results from the Supcr- 
Kamiokandc collaboration 111, as well as with values of 



e^r extracted from global fits to data from multiple ex- 
periments 



26 2 



In summary, this is the first direct search for non- 
standard interactions with high-purity samples of both 
neutrinos and antineutrinos. We conducted a simulta- 
neous fit to neutrino and antineutrino energy spectra of 
conventional Vp — > v T oscillations with an additional NSI 
matter effect. We found no evidence for non-standard 
neutrino interactions. 

This work was supported by the U.S. DOE; the United 
Kingdom STFC; the U.S. NSF; the State and Univer- 
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FIG. 2: The 68% and 90% allowed regions for the three pa- 
rameters in the combined oscillation and NSI flavor change 
model given in Eq. ((3|. To obtain slices for each parameter 
pair combination we marginalize over the third parameter as 
well as over the systematic penalty parameters. 
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